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Abstract 



We study constraints on type-II two Higgs doublet models at large tan (3 from 
LEP/SLD Z-pole data and from lepton universality violation in W decay. We 
perform a global fit and find that, in the context of Z decay, the LEP/SLD 
experimental values for lepton universality violation, R^, and A^ all somewhat 
disfavor the model. Contributions from the neutral Higgs sector can be used to 
constrain the scalar-pseudoscalar Higgs mass splittings. Contributions from 
the charged Higgs sector allow us to constrain the charged Higgs mass. For 
tan (3 = 100 we obtain the la classical (Bayesian) bounds of 

m H ± > 670 GeV (370 GeV) and 1 > m h o/m A o > 0.68 (0.64). 

The 2a bounds are weak. Currently, the Tevatron experimental limits on 
lepton universality violation in W decay provide no significant constraint on 
the Higgs sector. 

12.60.Fr, 12.15.Lk, 13.38.Be, 13.38.Dg 



Typeset using REVTpjX 

*electronic address: lebedev@quasar.phys.vt.edu 
^electronic address: loinaz@alumni.princeton.edu 
^electronic address: takeuchi@vt.edu 

1 



I. INTRODUCTION 



Perhaps the most important unanswered question in particle physics today is: "What is 
the nature of electroweak symmetry breaking?" The Standard Model (SM) incorporates the 
simplest mechanism: a Higgs sector consisting of a single self-interacting scalar SU(2) dou- 
blet of hyper charge Y = 1. Upon breaking of electroweak symmetry, the physical spectrum 
of the SM Higgs sector consists of one CP-even neutral Higgs particle. Current experimental 
data do not definitively contradict the SM, but persistent deviations in precision electroweak 
data from SM predictions on the edge of statistical significance tantalize with the possibility 
of new physics. This, together with various theoretical prejudices which suggest that the 
SM cannot be a complete theory, motivates the detailed study of alternative scenarios of 
eletroweak symmetry breaking (EWSB). 

The Two Higgs Doublet Model (2HDM) |IJ is the most straightforward extension of the 
EWSB mechanism of the SM. The theory proposes a pair of scalar SU(2) doublets, both with 
hypercharge Y — 1. Depending on the version of the 2HDM, these scalars may couple in 
various ways to the quarks and leptons. After electroweak symmetry is broken, the spectrum 
of the Higgs sector consists of five physical Higgs bosons: two neutral CP-even scalars {Hq 
and Ho), a neutral CP-odd scalar (Aq), and a pair of charged scalars (H^. These particles 
could be detected via direct production at colliders, but their effects may also be visible 
indirectly, through their contributions as intermediate states in decay processes. 

In this paper we consider the indirect signatures of the 2HDM in flavor-conserving W 
and Z-decays through its contribution to decay amplitudes via loop corrections. We consider 
only Type-II 2HDM models, in which the I3 = | fermions couple to one Higgs doublet and 
the J 3 = — I fermions couple to the other. We also focus on the large tan/3 regionQ, in which 
the Higgs couplings to the down-type quarks and the charged leptons are enhanced^. This 
can potentially lead to observable (or constrainable) flavor dependent corrections in Z and 
W decay, especially for the third generation {b and r). 

One-loop corrections to flavor-conserving Z decays in the 2HDM have been consid- 
ered previously in Refs. including as a possible explanation for the now-defunct 
anomaly'. The Z-pole runs at LEP and SLD are complete and essentially all of the data have 
been analyzed. The anomaly' has disappeared only to be replaced by the l A^ anomaly' 
thus, it is timely to revisit the model. We perform a global fit to all LEP/SLD Z-pole 
observables, and we examine the competing constraints from lepton universality, Rb, and 
At, on the charged and neutral sectors of the model. In addition, we study constraints on 
the model from lepton universality violation in W decays, which have not been previously 
considered. 



Perturbativity of the b and t Yukawa couplings requires 0.3 < tan/? < 120 (see, for example fj). 

2 This model is often studied embedded in the Minimal Supersymmetric extension of the Standard 
Model (MSSM) j|, although we do not consider it in this context here. 
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II. LEPTONIC W DECAYS 



In this section we calculate the constraints on the large-tan (3 2HDM from lepton univer- 
sality violation in W decays. We use the Feynman rules and conventions of Ref. ||10|| . Our 
notation for the scalar and tensor integrals is established in Ref. [11]. 

The leading (in tan/3) one-loop corrections to the decay W~ — > r~v T are shown in Fig. [I]. 
The corresponding contributions to the amplitude are Q 



i2 / m r tan/? N 



-z^=^(Q)f(p) T P^(g) 

(la/h°) : sin 2 a 2(7 24 (0, 0, Q 2 ; 0, m h o,m H ±) 
(la/H°) : cos 2 a 2(7 24 (0, 0, Q 2 ; 0, m H o,m H ±) 
(la/ A ) : 2(724(0, 0, Q 2 ; 0, m A o, m H ±) 

(lb/h°) : sin 2 a B x (Q;Q,m h o) 
(lb/H°) : cos 2 a Bi(0; 0, m H o) 

(lb/A°) : B 1 (0;0,m A o) 
(lc) : 2^(0; 0, m H± ) 



(2.1) 



with Q 2 = m 2 ^. The tree level amplitude is the expression in the square brackets. For the 
diagrams involving h° and H°, we have dropped terms subleading in tan/50. In the above 
we have made the large tan {3 approximations: 



cos(/3 — a) ~ sin a 
sin(/3 — a) ~ cos a 
sin « 1 
cos /3 « 0. 



(2.2) 



Combining the above corrections, with factors of 1/2 for the wave-function renormaliza- 
tion diagrams (lb) and (lc), leads to a shift in the Wtv t coupling given by: 



9 



g 2 ( m T t&n(3\ 2 
2 V m w j 

sin 2 a |c , 2 4 (0,m /l o,m H ±) + - 5i(0,m h o) + -Si(0,771h±] 
+ cos 2 a ^C 24 (0,m H o,m H ±) + ^B^O, m H o ) + ^i?i(0,m^±^ 
+ |(72 4 (0,m A o,m^±) + -^51(0,777^0) + i73i(0,m^±; 



(2.3) 



3 In computing the one-loop vertex corrections in both W and Z decays, it is a good approximation 
to neglect light fermion masses in loops. 

4 The subleading contributions of the h° and H° diagrams combine with the diagrams involving 
the Goldstone bosons to give finite results. 
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where we have suppressed the external momentum dependence of the integrals for notational 
simplicity. Similar shifts to the Wfiu^ and Weu e vertices exist but they are suppressed by 
factors of (m^/m T ) 2 and (m e /m T ) 2 so we neglect them. 

The complete expression for the finite combination of integrals seen in the curly brackets 
of Eq. (|IO|) , namely, 



C(Q 2 ; mi, m 2 ) = C 24 (0, 0, Q 2 ; 0, m u m 2 ) + - Bi(0; 0, mi) + - Bi(0; 0, m 2 ), (2.4) 

can be found in the appendix of Ref. [[Ti l. However, for our purposes it will suffice to expand 
it in powers of Q 2 = mf^: 



2 



2 2 2m 2 m\ m\ 
m 1 +m 2 5 2 l n — 9 

mf — TO2 Wlj 



(4tt) 2 12(m? - m 2 ,) 2 
+ ■•• (2.5) 



G{x) = l + \ (P^-)\nx. (2.6) 



where 



2 VI -x 

Observe that the function G(x) is negative semi-definite so that the leading term is non- 
negative for all masses mi and m%. This term dominates the subleading term unless the 
splitting between mi and m 2 is less than about m^/2. In the limit that the masses are 
degenerate, the leading term vanishes and the expansion reduces to 

C(m 2 y ;m,m) = -^^ + ... (2.7) 

In the limit mi — > oo (the full expression is symmetric in mi and m 2 ) the expression becomes 

,9 ^ 11/1 ml\ , 

C(m 2 w ; mi, m 2 ) =► - — - (l + - In -| J + ■ • ■ (2.8) 

Thus, Eq. (|2.3|) appears to lead to non-decoupling of heavy particles. That is, when m^o 
and m H ± are taken to be large independently, the amplitude does not vanish. However, in 
the general 2HDM the mass eigenvalues and mixing angle are related in such a way that if 
m^o — > oo while the couplings and the W mass are held fixed, the Higgs masses and mixing 
angle approach the limit [1|]: 

m H o ~ m H ± ~ m A o ; m h o /m A o — > ; cos(a — (3) — > 0. (2.9) 

In this limit the amplitude vanishes and decoupling is obtained. This decoupling behavior 
can be understood as follows: In the large tan/3 limit, sin a — > and the two Higgs doublets 
do not mix. Since large tan/? implies v± — > 0, electroweak symmetry is unbroken in the $i 
sector. The leading tan 2 (3 diagrams are then due to the W~&i*<§>\ vertex and corresponding 



4 



wave function renormalization diagrams (i.e. with $5 an d $r m the loop). The sum of these 
diagrams vanishes at p 2 = as a consequence of the Ward identity, and thus the heavy Higgs 
bosons decouple.^ 

Setting mjifi = m#o = rriH± = m and neglecting the m h o contribution, the shift in the 
coupling, Eq. ( |2.3|) , is small and positive: 
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g 2 I m T tan/3' 
2 \ m w / 
g 2 ( m r tan /3 N 
2 \ m w t 
' gm T tan/3\ 



x 2((m 2 v ;m,m) 



m 



\v 



(An) 2 18m 2 



\ 247rm J 

Away from this limit, the shift in the coupling is negative. 



(2.10) 



&9t 
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1 g 2 ( m T tan/3' 



(4tt) 



x 



m w 

/ 9 

m 



sm 2 a G('-^)+cos 2 a G(^)+G{'-^ 



m 



h° 



ni 



HO 



< 0. 



m. 



m 



AO 



(2.11) 



So, the model predicts a negative 5g T , except in the limit that the Higgs mass splittings are 
small (< rnty/2). The magnitude of the shift is maximal for an extreme non-decoupling 
case in which the charged Higgs is much heavier than the neutral Higgses. In this case it 
reduces to 



Sg T 

g ~ (4tt) 2 4 



1 g 2 ( m T tan/3' 



m w 



ml 



(2.12) 



if we assign a common mass, m , to the neutral Higgses. 

The current bound on lepton universality violation in leptonic W decays from the D0 



Collaboration is [12 



— = 1.004 ± 0.019 (stat.) ± 0.026 (syst.). 

9e 



(2.13) 



The central value of 5g T is positive, which is not allowed when the leading G(x) term 
dominates. However, this fact is inconclusive since the experimental error is large. Using 



5 The leading tan/3 contribution of the hP boson does not exhibit decoupling by itself : it is 
proportional to tan 2 /3 sin 2 a lnm^ — ► lnm^ since sina ~ — cos (5+ O(m 2 z /m 2 A0 ) in the decoupling 
limit ||]. As the result is independent of tan/3, subleading tan/3 diagrams must be included to 
obtain the decoupling behavior. 
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m T (mw) = 1-777 GeV and 2mw/g = v = 246 GeV, and adding systematic and statistical 
errors in quadrature, we obtain from Eqs. ( 2.12 ) and ( p. 13 ): 



Wo J 



100 

tan/3 



1.2 ±9.7 



(2.14) 



Since G(x) negative semi-definite and invariant under x <-» 1/x, at tan/5 = 100 this leads 
to a la bound of 



m 



ml 



> -8.5 



which translates to 



or < 1.3 x 10 . 

m#± mo 



For smaller tan (3 the bound is even weaker. 

Similarly, if we assume the limit of Eq. (|2.10|) , the best-fit value of the common mass is 
(in GeV) 



f 100GeV\ 2 



m 



100 

tan/5 



17 ± 132 



(2-15) 



At la and tan/5 = 100, this translates into 

m > 8 GeV 

so, the bound is extremely weak in this mass-degenerate limit as well. Thus, even for 
tan/5 = 100 the current data gives no significant la constraint on the Higgs masses. 



III. CONSTRAINTS FROM LEP/SLD OBSERVABLES 

In this section we perform a global analysis of LEP/SLD precision electroweak data in 
the context of the large tan (3 2HDM. We calculate the linearized shifts in the Zff couplings 
from SM predictions, fit these shifts to the data, and use the results of the fit to constrain 
model parameters. 



A. Corrections to the couplings 

As in the W decay case, large tan f3 enhances the coupling of the Higgs sector to charged 
leptons and down-type quarks, but even then one only needs to consider the third generation 
fermions. Below we list corrections to Z — ► bb, rf, v T v T . 

The leading tan (3 corrections to the Z — > b^b^ are shown in Figs. and |3|. The ampli- 
tudes of these diagrams are: 
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-Z^Q)b R (p) lfl b R (Q) 



x 



g 2 ( nib tan /3" 
4 y m w J [ cos6 w 
(2a + 2b/ h ) : sin 2 o: 2(7 24 (0,0,Q 2 ;0 
(2a + 2b/ H°) : cos 2 a 2(7 24 (0,0,g 2 ;0, 

(2c/h°) : a 6i sin 2 a {(d - 2)(7 24 (0, 0, Q 2 ; m^o, 0, 0) 

-Q 2 C 23 {0,0,Q 2 ;m h o,0,0)} 

(2c/ H°) : a 6i cos 2 a {(d - 2)(7 24 (0, 0, Q 2 ; m H o, 0, 0) 

-Q 2 (7 2 3(o,o,Q 2 ;mHO,o,o)} 

(2c/ A ) : a 6i {(d - 2)(7 24 (0, 0, Q 2 ; m A o, 0, 0) 

-Q 2 (7 23 (o,o,Q 2 ;m A o,o,o)} 

(2d + 2e/a°) : 2h bR sin 2 a ^(0; 0, m h o) 
2hb R cos 2 a -Bi(0; 0, m H o) 
2h bR B 1 {0;0,m AO ) 

-4a H +(7 24 (0, 0, Q 2 ; m t , m H ± , m H ±) 
2h tL {(d - 2)C 24 (0,0,Q 2 ;m H ±,m t ,m t ) 

-Q 2 C 23 (0,0,Q 2 ;m H ±,m t ,m t )} 

(3c) : -2h tR m 2 C (0,0,Q 2 ;m H ±,m t ,m t ) 
(3d + 3e) : 4a fefl I?i(0; m t , m#±) 



(2d + 2e/H° 
(2d + 2e/A°) 
(3a) 
(36) 



(3.1) 



where 



hj = hf - Qf sin 2 6 



w 



(3.2) 



and Q 2 = m 2 z . The tree-level amplitude is the expression in the square brackets times h bR . 
As before, we have dropped terms subleading in tan/3. Combining these corrections, with 
factors of 1/2 for the wave-function renormalizations, leads to a shift in the right-handed 
coupling of the b to the Z given by 



6h bR = 6h?„+6h: 



N 



b R 



where 



(3.3) 



Sh 



N 



g 2 ( nib tan/3" 
4 V m w , 



x 



sin 2 a|2(724(0,m /l o,m j4 o) + ^B x (^,m }fi ) + -.B^O, m A o) j 



+ cos 2 aj2(7 24 (0,m H o,m A o) + ^B 1 (0,m H u) + ^B 1 (0,m A o) 
+h bL sin 2 a|(d - 2)C 2 4{m h o, 0, 0) - mlC^m^o, 0, 0) + Bi(0,m h o) 
+h L cos 2 a\ (d - 2)C 2 4(m H o,0, 0) - m 2 z C 23 (m H o, 0, 0) + B^O^ho) 



1 



+h b A (d - 2)C M (m A o,0,0) - m 2 z C 23 (m H o,0,0) + B 1 (0,m A o] 
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c g 2 I m b tan/3' 



x 



2 y m w J 
-h H + 1 2C^(m t ,m H ±,m H ±) + B x (m u m H ±] 

+h t A {d-2)C 2 4{m H ±,m t ,m t ) - m 2 z C 2 3(m H ±, m t , m t ) + B^m^m^] 



-h tR m 2 Co(m H ±,m t ,m t ) 



(3.4) 



As in the W decay case, these expressions can be well approximated by their leading terms 
in an expansion in m 2 z as long as the mass splittings among the Higgses are not small. Using 
the formulae from the previous section and from the Appendix, we find 



1 g 2 f nib tan (3 



2 r 



(4tt) 



niw 



sin a G I 



m 



h° 



cos 2 a G 



ni 



m 



m 



+ 



1 g 2 f nib tan/?" 



mr 



(An) 2 4 V m w 



m 



(3.5) 



H±. 



where the function G(x) was defined in Eq. 

F( X) = ^ (l + 

See appendix for details. 

The diagrams which correct the decay Z - 
and |3| with the replacements b R <-> b L and 
diagrams are 



and 



1 — x 



(3.6) 



Ol&l is the same as those shown in Figs. |2| 
-> The amplitudes of the neutral Higgs 



J / m;,tan/? N 



4 y niw 
(2a + 2o/a°) 
(2a + 2b/ H°) 
(2c/h°) 

(2c/H°) 

(2c/A°) 

(2d + 2e/h°) 
(2d + 2e/H°) 
(2d + 2e/A°) 
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cos 6 W 



x 



sin 2 a 2(^24(0, 0, Q 2 ; 0, m h o, m^o) 
cos 2 a 26*24(0, 0, Q 2 ; 0, m H o,niAo) 

-h bR sin 2 a {(d - 2)^24(0, 0, Q 2 ; m h o, 0, 0) - Q 2 C 2 3(0, 0, Q 2 ; m h o, 0, 0)} 
-a bjJ cos 2 a {(d - 2)^(0, 0, Q 2 ; m H o, 0, 0) - Q 2 C 23 (0, 0, Q 2 ; m^o, 0, 0)} 
-h bR {(d - 2)634(0, 0, Q 2 ; m A o, 0, 0) - Q 2 C 23 (0, 0, Q 2 ; m A o, 0, 0)} 

—2hb L sin 2 a £?i(0; 0, m^o) 
—2hb L cos 2 a -Bi(0; 0, m#o) 
-2/ i6£J B 1 (0; 0,m A o) 



(3.7) 



with Q 2 = m|. The charged Higgs diagrams lead to corrections proportional to 

/ m t cot (3 \ 2 
V m w ) 
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and are suppressed compared to the neutral Higgs diagrams by a factor of (mt/m b tan 2 f3) 2 ~ 
(7.6/ tan/3) 4 so will be neglected. The shift in the left-handed coupling of the b to the Z is 
then 



5h bL = ShZ + Sh 



N 



(3.8) 



with 



5h» = 



g 2 



m b tan/?" 



x 



sin 2 o;< 2(724(0, m^o, m^o) + -5x(0, m^o) + -_Bi(0, m^t/ 



+ cos 2 a< 2(7 2 4(0, m^o, m^o) + —Bi(0,m H o) + -Bi(0, m^o) 



sin 2 a|(d - 2)C 2i (m h o, 0, 0) - m z C 2 3(m h o, 0, 0) + 1^(0, ra^o) 
-/i 6fl cos 2 aj(d - 2)C 24: (m H o, 0, 0) - m 2 z C 23 (m H o, 0, 0) + 5i(0,m#o) j 
- 2)C , 24 (m A o, 0, 0) - miC^m^, 0, 0) + ^(0, m A o) 





(3.9) 



Again, in the approximation Q 2 — > 0, we find 



1 g 2 / m&tan/? N 



(4tt) 2 8 



sin a Gl 



h° 



cos a G 



H° 



m 



-Sh. 



N 



AO, 



m 



AO/ J 



(3.10) 



So in this approximation, the shift in the left-handed coupling of the b quark due to neutral 
Higgses is equal in magnitude but opposite in sign to the shift in the right-handed coupling. 

To estimate the corrections to Z — > uu, cc, we note that the Higgs couplings to u and c 
quarks are suppressed either by tan/5 or by small d and s quark masses. Thus, we neglect 
these corrections. 

The corrections to the r couplings to the Z can be obtained from those of the b couplings 
by the simple substitutions 



which lead to 



-Sh 



N 

TL 



K. = 

<"*?, = < - <>■ 



m t ,m b 
ht L , h tR 
h bL , h bR 



1 g 2 ( m T tan/3 N 



0, m T 
h TL 7 h TR 



(3.11) 



(4tt) 2 8 



m w 



sin a G I 



m 2 



+ cos a G i 



'm|o\ 



a , 



m 



AO . 



Note that the charged Higgs contribution is zero since m t is replaced by m„ 
F(fflJHi) = F(0) = 0. 



(3.12) 
and 
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The decay Z — > i/ T i/ r is corrected by the diagrams shown in Fig. |[ The amplitude of 
these diagrams are 



g 2 ( m T tan/5' 
2 V m w t 



—i- 



cos 9 W 

(4a) : h TR {(d - 2)(7 2 4(0, 0, Q 2 ; m H ±, 0, 0) 

-Q 2 C 23 (0, 0,Q 2 ;m H ±, 0,0)} 

(46) : /i H+ 2C 24 (0, 0,Q 2 ;0 
(4c + 4d) : 2K L B 1 {0; 0,m H ±) 

with Q 2 = m|, resulting in a shift of the neutrino coupling by 



5h 



g 2 I m T tan (3 



I'L 



h H ±^2C 2 i{0,m H ±,m H ±) + Bi(0,m H ±] 



+h TR { (d - 2)C 24 (m H ± ,0,0)- m 2 z C 23 (m H ± ,0,0) + B x (0, m H ±] 



As a consequence of F(0) = we find 



5h° L = 0. 



(3.13) 



(3.14) 



(3.15) 



To summarize, we have found that the non-zero shifts in the fermion couplings in our 
approximation (Q 2 = 0) are: 



1 g 2 ( rrib tan/?' 



m w 



c 



(4tt) 2 

| 1 9 2 f m b tan{3 \ 2 F ( m 2 ^ 



m 



AO, 



m 



AO 



(An) 2 4 \ m w 



\ m H±) 



5h 



N 

TR 



-6h 



N 

TL 



m 



m 



I ) <• 



(3.16) 



Since G(x) is negative semi-definite, the shifts in the left-handed couplings of the b and the 
r due to the neutral Higgs sector are both always positive while the shifts in the right-handed 
couplings are always negative, and they are all proportional to the same linear combination 
of G-functions. Also, since —1 < F(x) < 0, the charged Higgs sector produces only a 
negative shift in hb R of magnitude at most 

1 g 2 ( rrib tan/3\ 2 



(4tt) 2 4 



m w 
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B. Fit to the data 



We have identified the relevant vertex corrections to Z decay in the large tan/3 2HDM. 
Using the LEP/SLD data to constrain their sizes will let us constrain the ratios 



m h° m H° j m t 



TTljifi tri^o m H ± 

All the neutral Higgs corrections are proportional to each other, so we will use Sh^ as the 
fit parameter. For the charged Higgs correction we will use Sh^ R . 

In addition to the proper vertex corrections, the 2HDM corrects Z decay through oblique 
corrections which can be expressed as corrections to the p-parameter and the effective value 
of sm 2 9w- Since we will consider only ratios of partial widths and asymmetry parameters 
in our fit, the p-parameter drops out from our analysis and we need only consider the shift 
in sin 2 6\y which we will denote Ss 2 .f\ We will not utilize 5s 2 to extract information on the 
2HDM because oblique corrections are generically sensitive to other sorts of new physics as 
well. 

The shifts to the Zff couplings in the large tan f3 2HDM can then be expressed as: 

6h UeL = 5h UliL = 5h UrL = 
5h eLR = 5h ilLR = 5s 2 

Sh TL = Ss 2 + 8h^ L 
6h TR = 5s 2 - 5h? L 

Skur o = Skr.r o = -'SS 2 



Shd LR = Sh SLR = ^5s 2 



3 

< 

3 



1 ~ 2 ( m b \ r.JV , r.C 



Sh bR = -5^-(-J)5h^ + 8h^ R (3.17) 



The dependence of various observables on ShP, Sh^ R , and 5s 2 can be calculated in a straight- 
forward manner. For example: 



5A e _ 4h eL h eR (h eR 5h eL - h eL 5h, 



Ae « - K R 

'a 

-OS 



Ah eL h eR {h eR h eL ) , 2 



h 4 - 

U e L ,l e R 

= -53.5 5s 2 

where the coefficient has been calculated assuming sin 2 6w = 0.2315. Similarly, 



6 Similar techniques for isolating oblique corrections into one or just a few phenomenological pa- 
rameters to extract constraints on proper vertex corrections were used in Refs. [0Jl^ 15|. 
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5A T 

SA FB (e) SA FB (n) 

A FB (e) A FB ([i) 
SA fb (t) 



-53.5 5s 2 + 3.96 5^4 



-1075s 



2 



-1075s 2 + 3.96 Shi 
A fb (t) 

X D ST p 

on* = orir = _ g4 5g2 _ 2 gg +o m§h c + o.307 Sa s 

R e Rfj, L R 

SR 

= -0.84 5s 2 + 5.07 5/^ + 0.184 Sh? + 0.307 5a s 

R T 

"■^ b _nioox„2 -in o xuN , n «wuC 



SRr 



0.182 5s 2 - 10.3 5/4; + 0.652 5/i^ 
-0.351 5s 2 + 2.89 5/^ - 0.184 8ffi R 
-54.1 5s 2 + 3.43 Sh? L - 1.73 5^ 



A FB (b) 

6 4^ = -58.75s 
^fb(c) 

_ n fid X„2 , o AO xuN 



A b 
S_Ac 
A n 



0.681 5s 2 + 3.43 8h? L - 1.73 5/i^ 
-5.19 5s 2 (3.18) 



for sin 2 6w = 0.2315, m T {m z ) = 1.777 GeV, and m h {m z ) = 2.77 GeV. 7 We have introduced 
the parameter 5a s to account for the deviation of a s (mz) from its nominal value which we 
chose to be 0.120Q: 

a s (m z ) = 0.120 + Sa s . 

We fit the expressions in Eq. (|3.18|) to the differences between the LEP/SLD measurements 
and SM predictions shown in Table |. The corresponding correlation matrices of the data are 
given in Tables || and [II. The SM predictions listed are for a SM Higgs mass of 300 GeV. 
Changing the SM Higgs mass has a negligible effect on all fit parameters except 5s 2 which, 
as discussed above, we do not utilize except as a fit parameter. 
The result of the fit was 

6h? L = -0.00021 ± 0.00029 
Sh? = 0.0049 ± 0.0060 
5s% = -0.00069 ± 0.00019 



7 Letting mt>(mz) change by 10% in either direction will only change the final central value of Sh^ L 
by about a fifth of a a. 

8 The value of a s (mz) from LEP is determined by fitting the SM to Ri, (£ = e, /i, r). Since we are 
considering extra corrections to the R^s, we must let a s (mz) float in our fit. 
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5a s = -0.0007 ± 0.0051 (3.19) 



with the correlation matrix for the fit parameters shown in Table [IV]. The quality of the fit 
was x 2 = 18.4/(18 — 4). The largest contributions to the x 2 come from Ap^{b) (3.5) and 
(2.5) which means that the 2HDM corrections do not improve the agreement between 
the theoretical and experimental values of these observables. 

In Figs. |5], |6], and ^ we show how different observables constrain the parameters Sh^ L , Sh^ R , 
and 5s 2 . Since Sh^ L is the only parameter which breaks lepton universality, it is most strongly 
constrained by the ratios Ri {t = e, /z, r). This is evident from Figs. |5] and |. This places a 
tight constraint on the size of the neutral Higgs correction to the b quark observables.^] The 
charged Higgs contribution, Shf R , must then fit all the heavy flavor observables, but due to 
the small experimental error on Rb, it is also constrained to be small. In Fig. [7], one sees 
that the overlap of the Alr and Afb(&) bands prefers a value of Shb R of about 0.04, far from 
the SM point at the origin.^] However, the Rb band does not allow this deviation, leading to 
the large x 2 ' s f° r ^4lr and Aw(b) mentioned above. Note also that the large tan/3 2HDM 
predicts $h^ R < so it cannot account for the 1 Ab anomaly' even if the constraint from Rb 
were absent. In fact, since the best fit value of Sh^ R is still small but positive, the 2HDM is 
slightly disfavored by the data. 

The same can be said of 8h^ L \ Since the experimental value of R T is smaller than those 



for R^ and R e , it is easy to see from Eq. ( 3.18Q that the data prefer a negative value of 5h: 



However, the large tan (3 2HDM predicts 5h^ > 



C. Constraints on model parameters 

Since 5h^ L > and 5h^ R < for the large tan/3 2HDM in our approximation, the best 
fit values of Sh^ L and Shf R consistent with the model are Sh^ L = Sh^ R = 0, corresponding to 
the Standard Model case. Thus, the large tan/3 2HDM does not mitigate the Ab problem, 
nor does it even improve agreement between the theoretical predictions for Rb and R T and 
the experimental data.0 

In order to extract the limits on the Higgs mass ratios from Eq. ( |3.19|) , we have per- 



formed both classical and Bayesian statistical analyses, in the latter assuming a uniform 

TL 



prior probability for the parameter regions 5h^ > and 



9 This was pointed out by Hisano et al. in Ref. ||. 

10 This is the L Ab anomaly' mentioned in the introduction. See, for instance, Refs. [|| and |J. 

11 Similar behavior in the context of the MSSM with R-parity violation was observed in Ref. [13]. 
There, the preferred values of the fit parameters were again the opposite sign of what the model 
predicted, and moreover, more than one to two a away from zero. This was a manifestation of the 
Ab anomaly. In the model considered here, the Ab anomaly is not as manifest in the fit results, 
since there are fewer new physics parameters. 
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> 5hg R > 



1 g 2 ( mbtan/3 



(4vr) 2 4 



m w 



(3.20) 



(Recall that 



1 g 2 / m&tan/T 



(4tt) 2 4 



and -1 < F(ar) < 0.) 

The corresponding {68%} and 



confidence limits on the fit parameters are: 



classical : Sh% R > { -0.0011 } [ -0.0071 ; 

5h? L < {0.00008} [0.00037] 

Bayesian : 5hf R > { -0.0021 } [-0.0050; 

Shi < {0.00011} [0.00025] 



(3.21) 



Using 2m w /g = 246 GeV and m h {m z ) = 2.77 GeV, the bounds on Sh^ translate into the 
bounds on F(x), where x = m 2 /m 2 H ±: 



classical : F(x) > 
Bayesian : F(x) > 



37 

tan/? / 

' 52 
tan p 




(3.22) 



For tan/5 < 94, the entire range of F(x) is contained in the classical 95% confidence region 
(since — 1 < F(x) < 0). It is therefore difficult to significantly bound the charged Higgs 
mass using this method unless tan/3 is quite large. Choosing tan/3 = 100 for definiteness, 
we translate the bounds on Fix) into bounds on m#± 



classical : m#± > 
Bayesian : m#± > 



{670 GeV} [40 GeV] 
{370 GeV} [120 GeV] 



(3.23) 



Bounds on m H ± are plotted as a function of tan f3 in Fig. |[ 

Bounds on the neutral Higgs sector masses from constraints on Sh^ L are more model- 
dependent than charged-Higgs bounds, since Sh^ L involves the masses of all of the neutral 
Higgses and the mixing angle a. Generally, requiring the magnitude of Sh^f to be small 
constrains the scalar-pseudoscalar mass splittings to be small. To give a concrete example 
(as in Refs. ||[7l), let us consider the limit a = (3 and m^o = m#o, 



Sh 



N 

TL 



1 9 



1 / m T tan/? N 



(47T) 



G\ 



m 



h° 



m 



AO, 



In this approximation the {68%} and 
semi-definite) are 



lower limits on G 



a" 



(3.24) 



(which is negative 
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classical : G\ 



Bayesian : G\ 



rn 



h° 



m 



A°. 



m 



h° 



> 



> <- 



m 



A , 



22 

tan (3 / 

' 25 
tan (3 



47 
tan/3 

' 39 ^ 
tan/? 



(3.25) 



With the choice tan (3 = 100 we find that our bounds on G translate into bounds on the 
h°-A° mass splitting (choosing the branch of solutions with m^o/m^o < 1): 



classical : 1 > 
Bayesian : 1 > 



m h o 
rriAo 

itia° 



> {0.68} [0.43; 

> {0.64} [0.5i; 



(3.26) 



More generally, bounds on m^o / m^o are plotted as a function of tan (3 in Fig. || In the 
limit sin a ~ 0, the result is the same, except that m^o/m^o replaces m^o/m^o in the above 
expressions. 



IV. SUMMARY AND CONCLUSIONS 

We have analyzed the implications of the large tan (3 2HDM for Z decays and for lepton 
universality violation in W decays. For Z decays we find that the generic predictions of the 
model do not improve agreement between theory and experiment. Further, the LEP/SLD 
experimental uncertainty in the measurement of lepton universality is sufficiently small to 
place significant constraints on mass splittings in the neutral Higgs sector for large tan f3. 
Constraints from the b decay parameters are sufficient to place bounds on the charged Higgs 
mass that are increasingly strong for increasing tan f3. For instance, for tan (3 = 100 we 
obtain the la classical (Bayesian) bounds of 

m H ± > 670 GeV (370 GeV) and 1 > — > 0.68 (0.64). (4.1) 

For W decays, the experimental central value from D0 slightly disfavors the generic 
prediction of the model, but experimental uncertainties are too large to usefully constrain 
the charged-neutral Higgs mass splittings. 
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APPENDIX: FEYNMAN INTEGRALS 



The integrals we use here are defined explicitly in JTTJ. In the approximation p 2 = 0, the 
one-loop diagrams which appear in this work are proportional to the following expressions: 




(d-2)C 24 (0,0, p 



(4tt) 



2 -m s ,m f ,m f s j -m 2 z C 2 3 (o,0,p 2 ;m s ,m f ,m f 



m. 



X - ln^J +f(x) 



oc 2C 2 4 [0,0, p ;mf,m s ,m t 



(4vr) 2 



2 A -- |n 7 " 9,I) 



or 



oc 2C 2 4 (0,0,p 2 ;0,m Sl ,m s ^j 



2(4tt)< 



3 m 2 In m 2 — m 2 In m 2 

A s 1 ^1 *2 °Z 

2 m£ - mj 2 




oc m 



Co (0, 



,p ;m s ,mf,m f 



(4tt)< 



[f(x)+g(x) 




oc 5i (0; m/, m s ) 



(4vr) 2 



m 



-IA e -ln-^ l- 9 (x) 



where 



/(*) 



1 



4(1 - a;) 2 

!, 1 
-liax + — -r 

2 4(1 -a;) 2 



x — 1 — 2 In a; 



-(1 -x)(l -3x) + 2x 2 Inx 



for x = m 2 /m 2 . Note that the function F(x) appearing in Eq. (|3.6| ) is defined as 

F(x) = f(x)+g(x). 
For x — > 1 (degenerate scalar and fermion masses), 

n/ \ 1 X — 1 



x — 1 



g(x) 



+ 



(A6) 



(A7) 
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For x — > (the decoupling limit of heavy scalar masses), 

1 1 
f(x) « -lnx + -H 

g( x ) Ki -l\nx-~ + --- ~ -/(x). (A8) 

The function defined in Eq. ( j2.6| ) is symmetric under x <-> For x — > 0, oo, 

In a; 



2 

For x — > 1 (the decoupling limit of degenerate scalar masses) 



! + •••. (A9) 



G(*)~-^(*-l) 2 + ---. (A10) 
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TABLES 



Observable 


Reference 


Measured Value 


r /T'i'i"i'T'i) r > i • , • 

/r 1 1 1 kiH r rediction 


Z lineshape variables 






mz 





91.1872 ± 0.0021 GeV 


input 


Tz 





2.4944 ± 0.0024 GeV 


unused 


rr° 
"had 





41.544 ± 0.037 nb 


unused 


Re 


@ 


20.803 ± 0.049 


20.739 


Rfi 





20.786 ± 0.033 


20.739 


R T 





20.764 ± 0.045 


20.786 


^FB(e) 





0.0145 ± 0.0024 


0.0152 







f\ m r"~r i r\ r\r\-\ o 

0.0167 ± 0.0013 


0.0152 


^fb(t) 





0.0188 ± 0.0017 


0.0152 


t polarization at LEP 











0.1483 ± 0.0051 


0.1423 


A T 





0.1424 ± 0.0044 


0.1424 


SLD left-ri 


2jht asymmetries 






A LR 





0.15108 ± 0.00218 


0.1423 


A e 





0.1558 ± 0.0064 


0.1423 







0.137 ± 0.016 


0.1424 


A T 





0.142 ± 0.016 


0.1424 


heavy quar 


-c flavor 






Rb 





0.21642 ± 0.00073 


0.21583 


R c 





0.1674 ± 0.0038 


0.1722 


A FB (b) 





0.0988 ± 0.0020 


0.0997 


A FB (c) 





0.0692 ± 0.0037 


0.0711 


A b 





0.911 ±0.025 


0.934 


A c 





0.630 ± 0.026 


0.666 



TABLE I. LEP/SLD observables and their Standard Model predictions, 
predictions were calculated using ZFITTER v. 6. 21 [18] with rrit = 174.3 GeV 
and a s (mz) = 0.120 as input. 



The Standard Model 
H, m H = 300 GeV, 



18 





T"i 

m z T z 


o 

^had 


K e 


till 


rt r 


A FB (e) 


A ( \ 

-4fb(a) 


A ( \ 

^fb(t) 


m z 


1.000 -0.008 


-0.050 


0.073 


0.001 


0.002 


-0.015 


0.046 


0.034 


F Z 


1.000 


—0.284 


-0.006 


0.008 


0.000 


-0.002 


0.002 


-0.003 


o 

°had 




1.000 


0.109 


0.137 


0.100 


0.008 


0.001 


0.007 


R e 






1.000 


0.070 


0.044 


-0.356 


0.023 


0.016 


Rn 








1.000 


0.072 


0.005 


0.006 


0.004 


R T 










1.000 


0.003 


-0.003 


0.010 


A FB (e) 












1.000 


-0.026 


-0.020 


A FB {fi) 














1.000 


0.045 


^fb(t) 
















1.000 



TABLE II. The correlation of the Z lineshape variables at LEP. 





Rb 


R c 


A FB (b) 


A FB (c) 


A b 


Ac 


Rb 


1.00 


-0.14 


-0.03 


0.01 


-0.03 


0.02 


R c 




1.00 


0.05 


-0.05 


0.02 


-0.02 


A FB (b) 






1.00 


0.09 


0.02 


0.00 


A FB (c) 








1.00 


-0.01 


0.03 


A b 










1.00 


0.15 


A c 












1.00 



TABLE III. The correlation of the heavy flavor variables from LEP/SLD. 









5s 2 


5a s 


Sh? L 


1.00 


0.62 


-0.12 


-0.30 






1.00 


-0.22 


-0.63 


5s 2 






1.00 


0.25 


5a s 








1.00 



TABLE IV. The correlation matrix of the fit parameters. 
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FIGURES 





(d) (e) 

FIG. 2. One-loop neutral Higgs corrections to Z — ► 6_r6r. Diagrams which correct Z — > 
can be obtained by the interchange <-> 6r. 
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FIG. 3. One-loop charged Higgs corrections to Z — > b^bR. Diagrams which correct Z — > 
can be obtained by the substitution 6r — > &£, — > £r. 
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22 



-0.004 -0.002 0.002 0.004 

FIG. 5. The la constraints on 8h^ L and Sh^ from various observables in the 5s 2 = 5a s = 
plane. The shaded contours represent the 68% and 90% confidence limits. 
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24 



> 



600 
500 
400 
300 
200 
100 


50 60 70 80 90 100 

tan/3 

FIG. 8. Lower bounds on the charged Higgs mass vs. tan (3. The dot-dashed and dotted lines 
correspond to the Bayesian 68 % and 95 % confidence levels, respectively. The dashed and solid 
lines correspond to the classical 68 % and 95 % confidence levels, respectively. 
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50 60 70 80 90 100 

tan/3 

FIG. 9. Lower bounds on the scalar-pseudoscalar mass ratio vs. tan/3. The dot-dashed and 
dotted lines correspond to the Bayesian 68 % and 95 % confidence levels, respectively. The dashed 
and solid lines correspond to the classical 68 % and 95 % confidence levels, respectively. 
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